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Abstract 

We study /-electron lattice models which are capable of supporting oc- 
tupolar, as well as dipolar and quadrupolar order. Analyzing the proper- 
ties of the Fg ground state quartet, we find that (lll)-type combinations of 
the Fs octupoles Tf rl = Tj^+T^f+Tf are the best candidates for octupolar 
order parameters. Octupolar ordering induces Fs-type quadrupoles as sec- 
ondary order parameter. Octupolar order is to some extent assisted, but 
in its basic nature unchanged, by allowing for the presence of quadrupolar 
interactions. In the absence of an external magnetic field, equivalent re- 
sults hold for antiferro-octupolar ordering on the fee lattice. In this sense, 
the choice of our model is motivated by the recent suggestion of octupolar 
ordering in Np02- 

The bulk of our paper is devoted to a study of the effect of an external 
magnetic field on ferro-octupolar ordering. We found that octupolar order 
survives up to a critical magnetic field if the field is lying in specific direc- 
tions, while for general field directions, the underlying symmetry of the 
model is destroyed and therefore the phase transition suppressed even in 
weak fields. Field-induced multipoles and field-induced couplings between 
various order parameters are discussed on the basis of a group theoretical 
analysis of the Helmholtz potential. We also studied the effect of octupo- 
lar ordering on the non-linear magnetic susceptibility which satisfies novel 
Ehrenfest-type relations at continuous octupolar transitions. 



1 Introduction 

The nature of the 25K phase transition of Np02 is a long-standing mystery. 
The developments up to 1999 are reviewed in £Q. Np02 is a semiconductor 
with well-localized 5/ 3 shells, thus in principle relatively easy to understand in 
terms of crystal field theory and superexchange interactions. However, it was 
concluded that the phase transition at 25K cannot be modeled by any combi- 
nation of dipolar and quadrupolar ordering phenomena. This leaves us with the 
possibility that the primary order parameter is one of the octupolar moments 
,2: . Recent experimental evidence is successfully interpreted by postulating that 
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the primary order parameters are octupoles, and that the transition is ac- 
companied by the induced order of T$ quadrupoles 00]. 

In comparison to well-studied quadrupolar phenomena, the physics of models 
supporting octupolar order is less explored [S] . Motivated by recent suggestions 
that certain actinide and rare earth compounds undergo octupolar ordering 
transitions 121130! we study a lattice of /-shells with T 8 quartet ground states, 
assuming symmetry-allowed octupole-octupole and quadrupole-quadrupole in- 
teractions. Our choice of model is motivated by certain features of NpC>2 (and 
to a lesser extent, by Cei-^La^Be) physics, but we could not claim that we are 
offering a model for NpC>2. In particular, we consider only spatially uniform 
solutions though it is known that NpC>2 has antiferro-octupolar order following 
the four-sublattice triple-g pattern characteristic of the fee lattice 0] . Neverthe- 
less, our study of hypothetical ferro-octupolar ordering, in addition to having 
interest on its own, offers some insight into the physics of actual NpC>2. Some 
of the features of NpC>2 follow from the point group symmetry of the local 
Hilbert space, and from the mere fact of having symmetry breaking by IVtype 
octupolar order; we are able to account for these. We are, however, missing 
the interesting consequences of antiferro-octupolar intersite interactions, and 
Q ^ (0,0,0) space group symmetry. We plan to return to these questions in a 
sequel to our present work. 

Analyzing the possible order parameters carried by the r 8 ground state set, 
we find that the best choice of the local order parameter is the Ill-type combina- 
tion of T 5 octupoles. In contrast to the case of CeB 6 (which can also be modeled 
as an array of T$ shells), the fourfold degeneracy is lifted in a single continuous 
transition. We find that the primary ordering of octupolar moments induces 
Ill-type quadrupoles even in the absence of quadrupole-quadrupole interac- 
tions. However, allowing for a non-vanishing IVtype quadrupole-quadrupole 
interaction assists the octupolar ordering. These features are common to our 
hypothetical ferro-octupolar ordering model, and the actual antiferro-octupolar 
ordering of NpC>2. This holds in the absence of external (magnetic or strain) 
fields. However, our main interest lies in the study of field effects. 

The effect of external magnetic fields on ferro-octupolar transitions is ana- 
lyzed in detail. Since octupolar ordering is a mechanism for the spontaneous 
breaking of time reversal invariance, and the application of a magnetic field 
removes this invariance, it might have been guessed that the ground state de- 
generacy is completely lifted in a magnetic field, and consequently no octupolar 
symmetry breaking is possible. Such is indeed the case for fields applied in 
general (non-symmetrical) directions: the phase transition is smeared out by 
arbitrarily weak fields. However, we found that for magnetic fields lying in cer- 
tain planes, or pointing along specific directions, there is a remaining ground 
state degeneracy which is removed by a continuous symmetry-breaking transi- 
tion (or a sequence of two transitions). The critical temperature is gradually 
suppressed by increasing the field, eventually vanishing at a critical magnetic 
field. In simple terms, the explanation is the following: Magnetic fields will, in 
general, induce octupolar moments as a higher order polarization effect. How- 
ever, for special field directions, the field is not able to induce at least one of 
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the Vc, octupole moments. Then this moment can be generated by intersite 
interaction only, and it will become non-zero below a critical temperature. 

We stress that this latter part of our work rests on the assumption of uniform 
order, and we have not attempted to to apply similar considerations to the 
antiferro-octupolar phase. The essential difficulty is that an external magnetic 
field will, in general, turn the octupole moments from their original direction, 
and therefore the symmetrical zero-field triple-<f structure is expected to undergo 
a complicated distortion which is difficult to analyze. 

There is a great variety of multipolar moments induced by the concerted 
action of pre-existing uniform octupolar order and an external magnetic field. 
The existence of such polarization effects can be deduced from general symmetry 
analysis (Section^) and particularly Sections 15. II and l5.3|) : herein we follow and 
extend the approach by |Uj . Illustrative examples are provided by mean 
field calculations. 



1.1 Short review of NpC>2 

The experimental results on NpC>2 are known from Refs. ^lElElE]- We quote 
only the findings which are pertinent to our model study. 

Actinide dioxides have the CaF2 crystal structure at room temperature. The 
sublattice of the metal ions is an fee lattice. Np 4+ ions have the configuration 
5/ 3 , the corresponding Hund's rule ground state set belongs to J — 9/2. The 
tenfold degenerate free-ion manifold is split by the cubic crystal field, yielding 

a \Ta) doublet and the Ig } and ri^ ) quartets . Neutron spectroscopy has 
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Ig ) is well separated from the first 



shown that the ground state quartet 

excited state (the other quartet) 

5/ 3 states have both non-Kramers and Kramers degeneracy; the latter might 
have been expected to be lifted by magnetic ordering. There is, in fact, an ap- 
parently continuous phase transition at To = 25K, which was first observed as 
a large A- anomaly in the heat capacity The linear susceptibility rises to 

a small cusp at To, and stays almost constant below |12l 12"]. The observations 
were, at first, ascribed to antiferromagnetic ordering. However, diffraction ex- 
periments failed to detect magnetic order and Mossbauer effect has put the 
upper bound of on the ordered moment. Thus we can exclude dipolar 

ordering, and have to consider the possibilities of multipolar order. 

A 1?8 ground state quartet can carry a variety of local order parameters: 
in addition to the dipoles, also 1?3 and T$ quadrupoles, and T2-, I^-, and 
F^-type octupoles jZj. The nature of intersite interactions decides which of 
these will actually undergo an ordering transition. Though magnetic ordering 
is the commonly expected outcome, there is no rule to guarantee that dipolar 
ordering is the leading instability. In fact, for light rare earths (and also for light 
actinides) quadrupolar ordering often pre-empts, or preceeds, dipolar order. 
However, the possibility that octupolar order comes first, is quite novel 

The possibility of an explanation in terms of quadrupolar ordering has been 
examined with great care |14l HI H|. The ordering pattern has to conform to 
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the observed fact that the T < T phase preserves cubic symmetry, but this 
requirement can be satisfied by the triple-<f order of r$ quadrupoles [3]. Indeed, 
resonant X-ray scattering finds long-range order of the electric quadrupole 
moments. However, the unquestionable appearance of quadrupolar order cannot 
be the whole story. Quadrupolar ordering alone could not resolve the Kramers 
degeneracy |15j . and thus there should remain low-T magnetic moments giving 
rise to a Curie susceptibility; this is contrary to observations. Furthermore, 
muon spin relaxation shows that local magnetic fields, with a pattern suggestive 
of magnetic octupoles, appear below To 

Though there are cases (as in CcBg) when the Kramers and non-Kramers 
degeneracies of the Tg quartet are lifted in separate phase transitions, this is 
not the case for NpG^: quadrupolar and octupolar moments appear simultane- 
ously This alone suffices to show that they must be coupled, i.e., the octupolar 
moments must also be of T$ symmetry Still, the question may be posed 

whether the transition is primarily octupolar or quadrupolar. We confirm the 
proposal made in and 0]: the primary order parameter is the octupole mo- 
ment, and its ordering induces quadrupolar moments of the same symmetry. 
NpC>2 has antiferro-octupolar order. However, systems with a uniform polariza- 
tion of octupolar moments may, in principle, exist and their study is logically 
the first step in studying the nature of octupolar order, and its coupling to 
other kinds of order, and to external fields. Incidentally, we find that there is 
a phenomenological similarity between our findings and the observed behavior 
of NpC>2 (and also Cei-^LaaBg), the reason being that for Tg Hilbert spaces, 
the T§ octupolar states are non-magnetic, and magnetic susceptibility arises 
only due to transitions between octupolar levels. The mechanism relies only 
on the existence of an octupolar effective field, and it is basically the same for 
ferro- and and antiferro-type alignments. Naturally, the details of polarization 
phenomena depend on the kind of octupolar long range order, and we hope to 
return to the case of antiferro-octupolar order in a future work. 

In the absence of an external magnetic field, the thermodynamics of ferro- 
and antiferro-octupolar ordering is quite similar, as is also the case for more 
straightforward kinds of ordering. Our H = self-consistency equation for the 
uniform polarization is the same as the one for either of the sublattice polariza- 
tions in the triple-(f structure. Also the results for non-vanishing quadrupolar 
interactions are transferable since the coupled quadrupolar order is of the same 
T5 symmetry, and therefore the four-sublattice structure remains the same. 
For either ferro-octupolar, or triple-g antiferro-octupolar order we find that i) 
octupole-octupole coupling alone is sufficient to select a unique ground state; ii) 
the ground state carries also quadrupolar moment; iii) allowing for an additional 
quadrupole-quadrupole interaction does not change the character of the low-T 
phase until the quadrupole-quadrupole interaction exceeds a threshold value; 
beyond that, a purely quadrupolar transition is followed by octupolar ordering. 
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2 Octupolar moments in the Tg quartet state 

Neutron diffraction measurements indicate that the ground state is one of the 
two Tg quartets. Since the same irrep occurs twice, symmetry alone cannot tell 
us the basis functions: their detailed form depends on the crystal field potential. 
Exploiting the fact that the crystal field splittings are large, we neglect the 
higher-lying levels, and describe the phase transition within the ground state 
set. Since thus the sequence and separation of levels is of little consequence, 
we have arbitrarily chosen one of the quartets obtained by assuming a purely 
fourth-order potential 7i C ryst = O4 + 5C| as the ground state. We believe 
that this assumption has no influence on the main features of our results |17| . 
The general form of the basis states is (numerical coefficients will be given in 
Appendix |A"|) : 




It is apparent that the Tg quartet is composed of two time-reversed pairs, 
thus it has twofold Kramers, and also twofold non-Kramers, degeneracy. The 
above choice of the basis emphasizes the presence of J z dipole and 0% = | (3 — 
J {J + 1)) quadrupole moments, but of course it is not unique. In fact, the 
decomposition 

r 8 ©r 8 = ri g ©r 4 „©r3 ff ©r 5ff ©r2 M ©r 4 „©r 5M (2) 

[where g (u) indicates invariance (change of sign) under time reversal] shows 
that the subspace carries 15 different kinds of moments: three dipolar T4, two 
quadrupolar T3, three quadrupolar T5, and seven kinds (T2, T4,and T5) of oc- 
tupolar moments 0. Within the subspace we can rotate so as to get 
non-vanishing expectation values of any of the 15 potential order parameters. 
Or in other words, if there is an effective field —71(11) (where TZ may be any of 
the fifteen components), then the fourfold ground state degeneracy will be at 
least partially lifted. 

In general, intersite interactions affect all 15 moments. Lacking a microscopic 
mechanism of Np-Np interactions in NpC>2, we have to assume some form of 
the interaction, and argue from the consequences. 

Restricting our attention to octupolar ordering, we have to choose between 
T2, T4, and IV The possibility of T2 octupolar ordering was first suggested 
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by Santini and Amoretti j^j. In a later work, this choice was discarded because 
there would be no coupling to quadrupoles [3j . We may also add that T 2 ordering 
would still leave us with a twofold degenerate ground state. We can also exclude 
the 1?4 type octupolar moments, because symmetry allows their mixing with the 
magnetic dipoles, and dipole moments are excluded by experiments. 

This leaves us with the possibility of the ordering of r 5 -type octupole mo- 
ments [7] 

Tf = \(XJl-Jfh) ( 3 ) 

where the bars on the angular momentum operators mean the symmetrized 
combinations of operators J Z J% = JzJxJx + JxJzJx + JxJxJz, etc. Acting as 
a field, either of T£, Ty 8 , or Tf splits the Tg quartet into two doublets. The 
doublets carry magnetic moments. Thus assuming, say, 7~f -type ordering would 
still leave us with a residual degeneracy which should be lifted by a separate 
magnetic (dipolar) phase transition. 




Figure 1: Left: The direction dependence of the magnitude of the octupolar 
moment T(§, 4>) (see equation (0J). Right: The spectrum of T(#, <f>) in the [110] 
plane (<f> — tt/4). 

However, we may choose a different orthogonal set of octupole operators 
as order parameters. Fig. ^ (right) shows the spectrum of octupoles ^Hj 

T(0, 4>) = sin (cos tfTg + sin #ZJ?) + cos d 7? (4) 

for <j> = 7r/4. It appears that the ground state is always a singlet except for the 
special points d = and d — it (i.e., T@ which we discussed above). Further- 
more, the overall width of the spectrum varies with ??, reaching its maximum at 
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i9 = arccos or equivalent positions. Thus within the r 8 subspace, the 

three-dimensional pseudovector of octupoles is "longest" in the (111) direc- 
tion, or in equivalent directions. At any cj)) the maximum of the absolute 
value of the eigenvalues is taken as the magnitude of the octupolar moment 
(note that the spectrum is symmetrical about 0). This is the quantity shown 
in Fig. (left) . We see sixteen maxima. However, it is easy to check that the 
plotrange for (t9, (f>) gives two points for every direction, therefore the number of 
maxima is only eight. The directions (111), (Til), (111), and (111) are equiv- 
alent by cubic symmetry, and for each direction, the octupolc moment may be 
of either sign (Fig. ^ (right)). 

Thinking of the ordering as caused by an octupole-octupole interaction with 
cubic symmetry 

n oc = j oc YJ??A + T L T L + T t T L) ( 5 ) 

it is plausible that it will occur in one of the (111) directions. We may think of 
it as the cubic crystal field giving rise to an octupolar single-ion anisotropy with 
the (111) directions as easy axes. Thus our candidates for order parameters are 

T' 3 - T' 3 4- T' 3 4- T (i 

~ 1 x i 1 y ' x z 

111 z V z 

T P = _ r P , r f3 _ r f3 

z iTi x ~ v z 

T 0_ = _ T I3 _ T P +T r & \ 

The four minima seen in Fig. dbelong to ±Tf n and ±Tf T (23- We have checked 
that the ground state of either of these operators carries a Ts-type quadrupolar 
moment, but no magnetic dipole moment. 

We note that these single- ion properties would be useful for modelling NpC>2 ; 
however, the interionic interactions may still be chosen as either ferro-octupolar 
or antiferro-octupolar. Let us furthermore point it out that the four equivalent 
states © are ideally suited for constructing a four-sublattice ground state for 
nearest-neighbor antiferro-octupolar coupling on an fee lattice. This would cor- 
respond to the experimentally motivated suggestion of triple- if order by Caci- 
uffo et al. The triple- if octupolar ordering can then induce the observed 
triple- if structure of the T§ quadrupoles 

Cm = O xy + O yz + O zx , (7) 

etc., as a secondary order parameter. In the absence of a magnetic field, our 
mean field results are formally valid for either the ferro- or antiferro-octupolar 
case, while as far as field effects are concerned, we stick definitely to the former 
case. We have to refrain from making detailed comments on NpC>2 until we 
completed work on the magnetic properties. 
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Denoting the ground state of Tin by 
from Appendix 151 



we quote the numerical values 







and 



= A = -15.683 
<fo|0m|0o) = 5 = 8.019 

(<f>0\J X \<t> ) = (4>o\Jy\fo) = (fo\Jz\(i>o)=0 

(4> \o xy \ 4>q) = (00 \o yz \ 4>o) = (4>o \&xz\ 4>o) ■ 



(8) 



(9) 



One may be wondering whether the Y§ quadrupoles by themselves would 
like a different orientation than the one forced upon them by the octupoles. 
This is not the case: a calculation shows that the length of the pseudovector 
(O xy , O yz , O zx ) is the same in all directions. There is no single- ion anisotropy 
for the T 5 quadrupoles; picking the (111) solution is exclusively the octupoles' 
doing. 



3 The octupolar— quadrupolar model 

We assume the presence of Ts-type quadrupolar and octupolar interactions 

H = Hoc + Hquad (10) 

where H oc was given in (J5J and analogously 

^quad — ^quad ^ ^ (&i,xy& ^j,xy ~t~ Oi.yzOj.yz ~t~ @i,zx@ j,zx) ■ (H) 

For the sake of simplicity, we assume J qU ad/</oc _ 0. 

For ferro-octupolar coupling (J oc > 0), we may assume uniform (111) order. 
The mean-field single-site Hamiltonian is of the form 

Hmf = -Tf n (Tf n ) - jO in (Om) (12) 

where we have chosen the octupolar effective field amplitude as the energy unit, 
and j — Jquad/'/oc- Henceforth we assume that all effects arising from lattice 
geometry, and the detailed form of the interactions are included in J oc , J qua dj 
and hence also in j. 

We note that formally the same mean field problem arises by assuming 
antiferro-octupolar interactions, and postulating four sublattices with local or- 
der parameters as defined in © and the analogous quadrupolar moments. 

The temperature dependence of the order parameters T = (Tf^) and q = 
(Cm) is obtained by the numerical solution of the self-consistency equations 
derived from diagonalizing (|12|) in the basis QJ. The overall behavior is simi- 
lar to that found in dipolar-quadrupolar models used in the description of Pr 
compounds |^E2]; however, now octupoles play the role of dipoles. 
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Figure 2: The mean field phase diagram of the zero-field quadrupolar-octupolar 
model l|10|l in the quadrupolar coupling-temperature plane (j = J qua ,d/Joc and 
t = Tkft/ J oc )). The dashed and continuous lines signify first and second order 
phase transitions, respectively. Observe the regime of first-order transitions 
bounded by two tricritical points (marked by black dots). 



The dimensionless free energy belonging to (1 1 21) is 

F=)f 2 + \jq 2 - t \n(2exp(-Bjq/t)cosh(AT/t) + 2 exp(Bjq/t)) (13) 

where A and B were introduced in |JSJ, and t = k-aT/ J oc is the dimensionless 
temperature. 

Octupolar order (T ^ 0) induces quadrupolar moment even in the absence 
of a quadrupolar coupling, as we can see from setting j = in dJ-/dq = 

exp(AT/t)+exp(-AT/t)-2 
9 exp(AT/t) + exp(-AT/t) + 2 ' 1 ' 

The T — > limit is expressed in ©. It states that by construction, the (111)- 
type octupolar eigenstates carry (lll)-type quadrupolar moments. The same 
state of affairs prevails as long as T ^ 0. In the "para" phase above the 
transition temperature, all moments vanish. 

The continuous phase transitions of the model (|10fl can be described by the 
Landau expansion of the free energy (|13[l 

r fl A 2 \„ fj B 2 j 2 \ 2 IBjA 2 2 A 4 4 



2 4f / V 2 2t 4 t 2 * 96t 3 



2 



BjA 4 , 5 4 j 4 4 5 3 j 3 A , , 

^—rqT A + — V? —7—(t^T 2 + ... (15) 

24i 4 y 12i 3 q 12t 4 y v 7 



9 



where Tq is the non-critical part of the free energy. 

Critical temperatures are defined by the change of sign in the coefficient 
of either of the quadratic terms. Upon lowering the temperature, for small j, 
mixed octupolar-quadrupolar, while for large j, pure quadrupolar order sets in 
first. At intermediate j, there is a regime of first order transitions (Fig.^J. 

We consider first the weak-j limit. The critical temperature is t oc = A 2 /2. 
At f oc , octupolar moment appears as the primary order parameter, but there is 
also induced quadrupolar order. Minimizing T with respect to q, we get 

RA 2 T 2 

(16) 



U(jB 2 - t) 



Thus terms of order q 2 and qT 2 are effectively of 0(T 4 ). Minimizing with 
respect to T, the critical behavior of the octupolar and quadrupolar moment 



6B (18) 



A 2 - 8B 2 j V 2 

is characteristic of the mean field solution for primary, and secondary, order 
parameters. 

The coefficient 

A 4 AB 2 j-t 
96t 3 B 2 J - t 

of the combined fourth order 0(T 4 ) term of T changes sign at t = 4B 2 j. 
Equating this with the critical temperature, we identify the coordinates of the 
lower tricritical point as jtri,i = A 2 /SB 2 = 0.48, and itri.i ~ 123. The critical 
temperature is constant for j < jtri,i> F° r j exceeding jtri.i the transition be- 
comes first order. The nature of the coupled orders does not change, but they 
become more stabilized, and the common transition sets in at higher tempera- 
tures (Fig. EJ) . However, the ground state moments remain independent of the 
coupling strengths: qr^a = B, and Tt^q = A. Representative temperature 
dependences of T and q are shown in Fig. [3] 

At large j, the first instability is associated with the change of sign of the 
coefficient of the q 2 term: pure quadrupolar order sets in at t qua d = B 2 j. This 
critical line meets the boundary of first-order transitions at the critical end point 
j cn d ~ 2.75, tend ~ 177 (Fig. For j > j C nd there are two phase transitions: 
the onset of pure quadrupolar order is followed by the emergence of mixed 
octupolar-quadrupolar order at i oc . The lower phase transition is of first order 
up to the second tricritical point Jtri,2 ~ 3.75, ttri,2 ~ 185. For j < jtri,2, the 
onset of octupolar order is reflected in a discontinuity of q (Fig. right). For 
j > jtri,2, both transitions are continuous. 

Deep inside the quadrupolar ordered phase, the development of the octupolar 
order is essentially unaffected by what the quadrupoles are doing, apart from a 
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Figure 3: Octupolar (T) and quadrupolar (q) order parameters as a func- 
tion of t = k B T/J oc for Jquad/^oc = (left), J qua d/./oc = 0.75 (center), and 
■/quad/ Joe = 3 - 5 ( ri ght). 



weak effect on the transition temperature t oc (note in Fig. |2 that t oc saturates 
to a constant). Though in this regime, we cannot use Landau expansion to 
determine q, we may assume that it is near its ground state value B, and use a 
low-order expansion in T to obtain in the large- j limit 

km t oc = km ^exPWAoc) = ^ ^ M% 

j— oo i^oo exp(Bqj/t oc ) + cxp(-Bqj/t oc ) 

In familiar pkase diagrams of dipolar-quadrupolar models, the mixed order 
would be completely suppressed at Jquad/Jdipoic — > oo (see, e.g., Fig. 3 of 
Ref. In contrast, we find the finite saturation value (|19|l as J qU ad/Joc — > 

oo. The peculiarity of the situation depicted in Fig. [5] is the endurance of 
octupolar order even with infinitely strong quadrupolar coupling. Tke reason, 
as we understood earlier, is tkat in tke r s subspace tke T 5 quadrupoles are 
completely isotropic, tkus tkey can accomodate a reorientation of tke basis 
states witkout sacrificing any of tkeir rigid order. 

Even confining our attention to uniform states, tke effects of an external 
magnetic field are variegated: It may gradually suppress octupolar ordering, 
witkout ckanging its ckaracter (111 direction); it may split tke transition in two 
(001 direction); it may ckange tke ckaracter of octupolar order but still facilitate 
a pkase transition (11c direction); or it may completely forbid octupolar ordering 
(non-symmetrical directions). We will understand tkis in detail in Sections 
0] and |5] A straightforward ckaracterization of field effects in terms of tke 
magnetization curve and its derivatives (tke susceptibilities) is possible in tke 
H || (111) case only. Tkis is tke subject of tke next Subsection. 

We empbasize tkat our entire analysis of magnetic field effects is confined 
to spatially uniform states, and does not cover tke cases of supercell ordering, 
suck as tke experimentally observable antiferro-octupolar order of NpC>2 . Tke 
basic difficulty is tkat tke inter-sublattice angles of various moments may get 
ckanged by tke field; tkis effect will be treated in a subsequent work. 



11 







2 



0' 



T 



H 



Figure 4: Left: The (^Tiix) = T octupolar order parameter as a function of the 
temperature for H = 0, 0.5, and 1.0 (H || (111), H in units of <?^b)- Right: the 
T — H phase diagram of the J oc = 0.02/cb, J qU ad = model for H pointing in 
the (111) direction. The transition is continuous all along the phase boundary. 

3.1 Non-linear susceptibility: the H || (111) case 

In certain symmetry directions such as (111), the magnetic field merely acts to 
suppress uniform octupolar ordering gradually (Fig. 0] left). It appears that 
the phase boundaries can be scaled onto a common curve by introducing the 
field-dependent transition temperature 



This bears some similarity to the field-induced suppression of antiferro- 
quadrupolar order in PrFe4Pi2 |22) . 

Our starting point is the mean-field-decoupled hamiltonian 



where the notations follow i|12[) . Jin = (J x + Jy + Jz), and in the Zeeman term 
H is the reduced magnetic field. 

Multipolar phase transitions, even when non-magnetic, tend to have a strong 
signature in the non-linear magnetic response. The case of quadrupolar transi- 
tions has been extensively studied |2.'ij . To obtain analogous results, we expand 
the free energy corresponding to (|21|l 



T oc (ff) 



l-a H -H 2 -b H -H A ... 



(20) 



T oc (H = 0) 




(21) 





+ 2 cxp (Bjq/t) cosh (y H H/t)] . 



(22) 
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Here gn and yn are the two parameters of the Zeeman splitting scheme of the 
T 8 subspace (Appendix . The overall shape of the phase boundary in the 
t-H plane is obtained by expanding F(T, q, H) in powers of T, and identifying 
the coefficient of the T 2 -term 

r , tt t s _ 1 A2 sinh (g H H/t) 

' 2V ' ' 2 2g H H cosh (g H H/t) + cosh (y H H/t) ' { ' 

Solving C2(H,t) = gives a line of continuous transitions in the t-H plane 
(Fig. E| right). 

It is interesting that the octupole ordered phase can be suppressed gradually 
by a magnetic field (Fig 0] left). Dipolc ordering and octupole ordering are 
two independent ways to break time reversal invariance. However, octupolar 
moments are due to currents with zero total circulation, while dipole moments 
arise from non-zero integrated circulation. In a finite field H || (111), the 5/ 3 
ion must be able to sustain both kinds of currents simultaneously. 

•FfT, q, H) has to be expanded to 0(H 4 ) in order to derive both the sus- 
ceptibility Xi and the non- linear magnetic susceptibility \3- We do not give the 
detailed formulas here, but discuss the terms giving rise a quadratic shift of the 
transition temperature in (121)1 

T{T,q,H) » F(T,q,H = 0) 

+ g -^ L B j qH* - A.H*A*Bj0* 

9 2 H + 3y 2 H A 2 T 2 H 2 (24) 

48i 3 V ' 

The first term in the second line describes field- induced T$ quadrupoles. The 
general nature of octupoles would allow the presence of a TH term; it is the 
peculiarity of the (111) direction that it does not appear here. We have omitted 
field-induced terms of the non-critical part of the free energy; they have to be 
included when calculating the susceptibilities. 

Further calculation is analogous to that given for the H = case. Minimizing 
with respect to q gives for the secondary order parameter 

gw - iF^ '( A ' 7a + ^-^- (25) 

Replacing this back into (|24|) . we can determine the optimum value of the pri- 
mary order parameter T. Here we quote only the result for the quadratic shift 
of the transition temperature (cf. Equation l|2()|l) 

aH = 6^(^-2fflj) [(A2 + ■ (26) 

One of the contributions to an vanishes at the tricritical point (j — -> A 2 /8B 2 ), 
while the other remains finite. 
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Figure 5: Linear susceptibility (left), temperature derivative of linear suscepti- 
bility (center) and nonlinear susceptibility (right) as a function of temperature 
for magnetic field parallel to (111) direction (J oc = 0.02fcB, ^quad = 0). 

Representative results for the linear susceptibility x = —d 2 J 7 /dH 2 , and the 
third-order susceptibility \3 = —d^F/dH 4 , are shown in Fig. El The octupolar 
transition appears as a cusp in x (Fig. El left). The cusp can be also represented 
as the discontinuity oidx/dT (Fig. El middle). The non-linear susceptibility has 
a discontinuity from positive to negative values (Fig. El right). These anoma- 
lies are related to each other, and the specific heat discontinuity AC, via the 
Ehrenfest-type equation 24 



The derivation [25] of J23) relies only on the fulfillment of to order H 3 . 
In particular, bn does not come into (|27H . The relationship (|20|l is often found 
for the critical temperature of transitions to non-ferromagnetic phases like anti- 
ferromagnets, spin-gapped phases, quadrupolar order, etc. Octupolar ordering 
belongs to this class of transitions. 

Though the example shown in Fig. El was for J qua d = 0, the relationship 
<|27l) holds everywhere along the lines of continuous phase transitions shown 
in Fig. |21 As long as we are dealing with ordinary second order transitions, 
Landau theory would be consistent with all the discontinuities appearing in 
<|27|) being finite. However, approaching a tricritical point AC — > oo, and l|27|l 
allows several scenarios. We note from <|26ll that generically yn ^ 0, thus ajy 
remains finite, and then the simplest expectation is that the divergence of AC 
is matched by that of Adx/dT. Such is indeed the finding for our standard Tg 
subspace specified in Appendix^ The same holds for T 8 subspaces derived from 
a combination of fourth-order and sixth-order crystal field potentials. However, 
at one particular value of the ratio of the sixth-order and fourth-order terms, the 
Zeeman spectrum consists of a doublet and two singlets, i.e., the yn = case 
is realized. For this special model 123 a# — » as one approaches the tricritical 
point, and the field dependence of the transition temperature is purely quartic 
<|20|) . Here, a peculiar form of the Ehrenfest relation can be derived, an — > 
cancels the mean-field divergence of AC, and at the same time A%3 — > 0. The 




(27) 



14 



discontinuity of (px/dT) is now balanced by that of the fifth-order non-linear 
susceptibility 

^A X5 = & ka(|). (28) 
It should be interesting to find a situation where l|28|) is experimentally testable 

EE] 

4 Suppression of ferro-octupolar order by mag- 
netic field 

Next, we consider the effect of a finite magnetic field of arbitrary orientation 
on a system of interacting T 5 octupoles. Henceforth, our mean field arguments 
will be based on a simplified version of l|21|) 

T~t = Hoc + Ti-z = Hoc — H-J (29) 

where TC oc is taken from JSJ, and </ qU ad = 0. Here, as in Section [SJ we will 
confine our attention to uniform states. An essential extension of our argument 
would be needed to cover the case of Np0 2 . 

First, let us reconsider the case of a field H || (111) (Fig. The nature of 
the octupolar order parameter is not influenced by the field, only its saturation 
value, and the transition temperature, are scaled down. Conversely: since there 
is a T — phase transition at a critical field H CI (Fig. |SJ left), there must 
exist (in mean field theory) a finite-T ordered phase at H < H CI (Fig. 0] left). 
Generally, the nature of the field direction dependence of octupolar ordering 
can be studied by confining our attention to the ground state (T — 0). First, 
we use mean field theory; later, we give general symmetry arguments. 

Let Eq be the minimal eigenvalue of the mean field hamiltonian which con- 
tains both the external magnetic field, and the octupolar effective field 

£o«T H » = (H) = - ($o |7h| $ ) < 3k) - H- (J) (30) 

where |$o) is the interacting ground state, and () denotes expectation values 
taken with l^o)- The energy unit is like in Eqn. I|21|) . We have to minimize 

E({T H )) = ±{T u ) 2 + E a ((T Ii )) (31) 

with respect to (3k). In general, 7k is not pointing in the same direction in the 
T$ space as the zero-field T = (Tf^), but neither is it collinear with H; it has to 
be chosen in an optimization procedure, observing the symmetry lowering due 
to the magnetic field. 

Fig. H3 shows results obtained by minimizing E(T) with respect to T = 
(3m)) our original choice of order parameter. For H || (111), the field does not 
introduce any inequivalence of x, y, and z, thus the above choice of the order 
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Figure 6: Left: Field-induced ground state transition from octupolar order to 
the disordered state for H || (111). Right: Expectation value of the octupole 
moment in the ground state as a function of magnetic field for H || (123). For 
fields pointing in non-symmetric directions, there is no sharp phase transition. 

parameter is optimal, and the second order transition seen in Fig. (left) is 
genuine. 

Fig. (right) shows the self-consistent solution for T for fields H || (123). 
(123) is taken to represent general non-symmetric directions. We find behavior 
characteristic of smeared-out phase transitions (the marked upward curvature 
at H ~ 1.5 shows where the phase transition might have been; clearly, intersite 
interactions are important for H < 1.5, while their effect is negligible in the 
high-field tail). 
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Figure 7: The Landau-type ground state energy expression has symmetric or 
asymmetric minima depending on whether the field is applied in a symmetry 
direction (H || (111), upper curve), or non-symmetry direction (H || (123), lower 
curve) . 

The reason for the discrepancy between the H || (111) and H || (123) behav- 
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ior becomes clear from plotting the Landau-type ground state energy density for 
different field directions (Fig. [TJ). For H || (111), equivalent minima remain at 
the positions ±7o, thus the system can pick one of these in a symmetry breaking 
transition (upper curve). On the other hand, for H || (123), the two minima are 
not equivalent, the ground state remains always on the right-hand side. There is 
no symmetry breaking transition though the T-dependence may be non-trivial, 
showing the shadow of the phase transition which might have happened. 

One might object that for H || (123), T x , T y , and T z are no longer equivalent, 
thus the optimal mean field solution should be sought in the form 



This is true but we do not have to make the considerable effort of a three- 
parameter optimization. We will bring general arguments to show that the 
solution would be like that in Fig. (right), whatever 7h is chosen. It remains 
true that for general (non-symmetric) field directions, the ground state of l|29[) 
is non-degenerate and therefore no symmetry breaking transition (in particular, 
no continuous octupolar ordering transition) is possible. 

For some symmetry directions which are not equivalent to (111) (e.g., for 
(001)), the character of the solution is different from either of those shown in 
Fig. |H| We will discuss these later. 



5 Symmetry analysis of field- induced multipoles 
5.1 The high-field limit 



We learn from Fig. that at sufficiently high fields the ground state is deter- 
mined by the external field only: either because ordering has been suppressed, 
or because there was no transition to begin with. The following analysis of the 
field-induced multipoles does not rely on the mean field approximation, but each 
of the cases will be illustrated by a mean field calculation. 

The octupole operators are third-order polynomials of J x , J y and J z . The 
T5 octupoles can be expressed in terms of dipole and quadrupole operators [H] 



7k = r x % + r y T y + r z T z . 



(32) 



t£ = {\ol + \ol)j x + ^{o zx j z -o xy j v ) 

Ty = {- l -Ol + \ol)J y + 2 -{O xy J x -O yz J z ) 

1 2 
T z = ~ 2^2"^z "^i^yz'^y @zxJx) 



(33) 



where the quadrupoles are well-known quadratic expressions 




O. 




1 
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O yz — ^ {JyJz + JzJy) 

O zx = ^(J Z J X + J X J Z ). (34) 

Q33J) contains exact identities respecting the non-commutative nature of the 
operators. However, Tf etc. are themselves defined as symmetrized expressions 
(|3J), so it must be true that the order of the operators on the right-hand side 
cannot really matter. In fact, there is an arbitrariness in the representation 
: it would be also true that 

Tf = -ppz (35) 

or 

= I (Pv~^v~ ~ OZX) • (36) 

Similar relationships can be listed for the first two lines of <|33|l . This suggests 
that the relationships l|33|) can also be interpreted in terms of c- numbers, i.e., 
classical polarization densities Such considerations are valid for either 

field-induced, or interaction-induced multipole densities. This enables us to use 
relationships like (|33() in Landau expansions. 

First, we discuss field-induced densities. The basic idea is that an exter- 
nal field induces (J) = ((J x ), ( Jy), ( Jz)) | H , and this gives rise to induced 
quadrupoles (O xy ) = (J x )(J y ) as a second-order effect, and induced octupoles 
as a third-order effect, etc. If an octupole component is field-induced, it can no 
longer play the role of the order parameter of a symmetry breaking transition. 
The question is, can it happen that certain octupole moments are not induced 
by the field. 

l|3"5|l and l|3"6"| are still separate operator identities, but they must have the 
same classical meaning when J x , etc. are treated as c- numbers. Indeed from 



T/ = O yZ J y -0 ZX J X ^((Jy)(J Z ))(J y )-((J Z )(J X ))(J X ) 

= ({J y ) 2 -{J x ) 2 ){J z )^(Hl-Hl)H z (37) 
which is the same that we would have obtained from (|35|l . Similarly, 

t£ = a^7T-o^;^((j,)<j :c ))<j 2 )-((j.)<j y ))<j y ) 

= {{J z ) 2 -{J v ) 2 ){J x )<x{Hl-Hl)H x , (38) 



and 



Zf = O xy J x -O yz J z <x{H 2 x -H 2 z )H y . (39) 



We can also argue in the following manner. Higher-order polarizations in a 
magnetic field give rise to the following general 77-dependence of the energy 

£{H) ~ £{H = 0)-^H 2 - ^H 4 ... (40) 
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The lowest order time reversal invariant expression containing T is TH, thus 
the coupling of octupolar moments to fields may appear in terms from 0{H A ) 
upwards. If it does, the minimal eigenvalue of the mean field energy l|3(J|) will not 
be symmetrical under the sign change of octupole moments: Eq(7~) ^ Eq(—T~). 
Non- equivalent minima like in Fig. [7J (lower curve) mean that there is no sym- 
metry to break, a phase transition is not possible. However, for fields in the 
special directions discussed above, there is no field-induced T 5 octupole |3"U) . 
the ±T minima of Eq(T~) remain equivalent (Fig. [7] upper curve) and sponta- 
neous symmetry breaking remains possible. The eventual merging of the two 
minima is no longer a question of symmetry, but of field intensity; a sufficiently 
strong field will suppress octupolar (or any other) order, and produce a unique 
polarized state for any field direction (Fig. |SJ right). 

In what follows, we calculate the induced T 5 octupoles using lj3Tjl - il3T)|) for 
several field directions, and discuss the possibility of symmetry breaking tran- 
sitions. 

5.1.1 Non-symmetric directions 

First let us observe that a field pointing in a general direction will give non- 
zero values for 7~f , T@, and T@ . Since the field induces all T 5 octupoles, there 
remains no degeneracy to be lifted, no symmetry breaking transition is possible 
pH] . This corresponds to the situation in the right-hand panel of Fig. 

5.1.2 H || (111) 

Taking now H || (111), we find T£ = = T/ = 0, and so also Tf u = 0. 
The field does not induce T5 octupoles, and therefore a symmetry breaking 
transition is possible. Furthermore, as we remarked earlier, the x, y, and z axes 
play equivalent roles, and therefore the choice of the order parameter T = (T{ xl ) 
is correct. The situation corresponds to Fig. H3 (left). 

5.1.3 H || (001) 

Next consider H || (001). Also here we find Tj? = = Tj> = from Eqs. (|57|l - 
and therefore the possibility of continuous phase transitions. However, it 
is intuitively clear that the z-axis is inequivalent to x and y, and therefore the 
order parameter may be either T z , or some linear combination of T x and T y . 
We have to perform a two-parameter minimization using the suitably modified 
form of (Jjnj 

E ({T Z ),{T ± }) = -J ((T x + T y ){T ± ) +%(%)) (41) 

with (T x ) = (T y ) = (71). Like in Fig.[7J we expect that the ground state energy 
functional has degenerate local minima: at low fields, we find a pair of these 
as a function of (T z ), and another pair along the (T x -\-T y ) direction (the latter 
choice is arbitrary in the sense that we could also have taken (T x — Ty)) (Fig. [SI 
top left). These two pairs of minima are not symmetry-related, as it is also 
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Figure 8: The contour plot of the ground state energy functional in the 
(T/ + T/)- {%) plane for H = (0,0, H) magnetic fields H = 0.2 (top left), 
H = 0.42 (top right) and H = 0.8 (bottom). 

shown by the fact that at intermediate fields, only the {T x + T y ) ^ minima 
survive (Fig. [51 top right). At high fields, the ground state is non-degenerate 
with (% + T y ) = (T z ) = (Fig. El bottom). 

The corresponding sequence of two second-order ground state transitions is 
shown in Fig. EJ (left) . (T£ +T y 8 ) > (or alternatively, (T£ -T/) > 0) order 
develops at higher critical field H>, with (Tf ) = 0. Upon reducing the field to 
a lower critical value H<, a second symmetry breaking occurs. In the low-field 
phase H < H<, {T£ + T?) ? and also (Tf) ? 0. (T/) ? (Tj? + T/)/2 as 
long as H > 0; the T\\\ order ((Zjf) — (T@) — (Tf)) appears continuously as 

5.1.4 H || (11c) 

H || (11c) (c f 0) induces (TJ? -Tf) + 0, leaving (T/ + 1$) = and (Tf) = 0. 
There is a remaining octupolar degeneracy which is is lifted in a single contin- 
uous transition, where (T/) ^ 0, and (TJ? + Ty 3 ) ^ appear simultaneously 
(Fig.El right). For c = 0, (i.e., H || (110)) Tf , and +T y 3 can order indepen- 
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Figure 9: Left: The field dependence of (Tj> + T p )/2 and {%) in a field H j| 
(001). Right: Octupolar components for H || (11c), c ^ 0. 



dently, like in the case of H || (001). 

5.1.5 Field direction dependence: Summary 

We have discussed field directions which do not subtend a too large angle with 
(111), thus it holds that the limit H — > picks the (111) ground state. For 
other field directions, the limit H — > may give one of the other ground states, 
e.g. (Ill) type order (see (©)• 

The previously discussed special directions which allowed a symmetry break- 
ing transition, were all lying in the plane with normal vector ft = (1, — 1, 0). Be- 
cause of the cubic symmetry, the behaviour is the same for magnetic fields lying 
in planes with normal vectors ft = (1,1,0), ~rt — (1,0,-1), ~rt = (1, 0, l),~ft = 
(0, 1, —1), and ft = (0, 1, 1), only the ordering phases change correspondingly. 
These six planes intersect along the directions [111], [Til], [ill] and [ill] 
fFig. I1U[1 . Any direction outside these planes excludes the possibility of a con- 
tinuous octupolar transition. 

5.2 Field-induced multipoles 

We may also regard the problem of field-induced octupoles as a special case 
of field- induced multipoles in general 0|B]. It is best to begin with quadrupoles. 
When the magnetic field is zero, the Ts-type quadrupolar moments O xy , O zx , 
O yz arc equivalent by cubic symmetry. Switching on an external magnetic field 
with general direction Tl\\Jt, where it = (k, A, /i), the quadrupolar operator 
along the ~rt direction can be given as 

Q(~rt) ~ 3(rf J ) 2 — J(J + 1) = 3kX(J x J v + J y J x ) + 3k/j(J x J z + J Z J X ) 
+3\n(J x J y + J y J x ) + n 2 Jl + X 2 Jy + ^Jl) - J(J + 1) 
= 3k\O xv + 3KfiO zx + 3AfiO yz + (/i 2 - A 2 ) (J 2 - J 2 ) 
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Figure 10: Field directions lying in any of the planes shown allow a continuous 
octupolar ordering transition. Special rules hold for the lines of intersection, 
and other high-symmetry directions. 



+ ( K 2 - f)(Jl Jl) + (A 2 - , 
+ (k 2 + A 2 +^ 2 -1)J(J+1). 



(42) 



This means that the O xy , O zx and O yz r 5 -typc quadrupolar operators are 
no longer equivalent, the quadrupolar moment is distorted along the external 
magnetic field direction. Above, we obtained the components of the quadrupolar 
operator in the rotated new basis. The T§ quadrupole moments O xy , O zx , O yz 
are proportional to n\ ~ H x H y , k\i ~ H X H Z , X^i ~ H y H z . The next three 
quadrupolar terms, which are linear combinations of the two r,3 quadrupoles, 
are shown because they appear in the expressions of the T 13 octupoles: J z —J y ~ 



H^-Hy, J x -J z ~ H% — H%, Jy-J x ~ Hy — H%. The last term is an invariant 



number. 

5.3 Ordering in external magnetic field 

Now we consider arbitrary field intensities. At sufficiently low fields, symmetry 
breaking transitions are possible. However, the magnetic field lowers the sym- 
metry of the system in a peculiar way, and gives rise to couplings between order 
parameters which would be independent in the absence of a field. The nature 
of these couplings depends sensitively on field direction. 

Our previous discussion was about the ground state energy Eo(y, H), where 
H is the external magnetic field, and y stands for all other variables. The 
magnetic moment, or in our case J, is obtained as J = — (dEo/dH) y . 

In what follows, we prefer to use the Hclmholtz free energy Q which is related 
to .Eo by the Legendre transformation Q = E n + J-H. The magnetic field will 
be expressed as 




(43) 
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The generalized Hclmholtz free energy can be expanded in terms of the 
components of the symmetry-allowed multipolcs 7 

S G{Jx 7 Jy i Jz i @2 i ®2 ) ®xy 7 ®yz j O zx 7 Txyzi T x ,Ty ,T Z ,T y ,T^) 

= £ r^aiy..) (44) 

We have to go over the list of all possible product representations spanned 
by the order parameter components, and identify the bases for the identity 
representation Ti; these are the invariants 2(Ti<g)Tj...) (assuming that they are 
also time- reversal invariant). 

It is obvious that the present argument is valid only for phases with uniform 
order, q ^ Fourier components of the multipole densities should be included 
in a Landau theory of modulated states, such as the antiferro-octupolar phase 
observed in NpC>2. 

Returning to l|44|) : to construct a Landau theory, we would need all the 
invariants up to some specified order, but to derive H, it is enough to consider 
those which contain the components of J. Their general form is T^T^Tj...) — 
JV, where the components of V give the basis of r 4tl (g and u refer to parity 
under time reversal). 

We arrange the invariants according to the number of factors in the un- 
derlying product representation. For the present purposes, we will call this 
number the order of the invariant |32| . The second order invariants containing 
J are JJ and J-T a . Third order invariants arise from TztCgTsigirs, r^raCgTs, 
r4(g)r 2 £g>r5, and T4(g)r4(g)r5. To take the simplest example, consider r4(g>r2(g>r5, 
for which 

z(r 4 ®r 2 ®r 5 ) = J x o yz T xyz + J y o zx T xyz + J z o xy T xyz , (45) 

and the corresponding term of V is 

V(r 2 ®r 5 ) = {O yz T xyzl O zx T xyzl O xy T xyz ) . (46) 

Taking e.g., the z-component, we find that the magnetic field couples to O xy %c yz - 
One possible interpretation is that, in the presence of H || (001), O xy -type 
quadrupolar moment induces the octupole T xyz 0. Alternatively, ^j^-type 
octupole order would induce O xy quadrupoles. 

In our further discussion of third order invariants, we confine our attention 
to those which have a bearing on the question of octupolar order, i.e., one of 
the factors is rs u . As for r4(g)F4(g)r5 invariants, since one of the T4 has to give 
J, i.e., it is r4„, the remaining T4 must be F4 9 . The lowest order T4 9 multipole 
is a hexadecapole. However, within our Ts subspace, hexadecapoles cannot be 
independent of the first 15 multipoles, thus the formally third order expression 
would have to be rewritten as a fourth-order invariant. We generally neglect 
terms of fourth order, and seek to draw conclusions from the genuinely third 
order terms. These belong to T4 U (E)T 5u (E)r 5g which gives 

i(r± u ®r 5u ®r 5g ) = u-OxyTy + o zx rf) + j y (-o yz Tf + o xy r£) 

+J z (-O zx T£ + O yz T y 3 ), (47) 
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and to r4„<g)r5 U (g>r 3£ , which gives 

i(r 4u ®r 5u ®r 3s ) = - 1 - j x (o° 2 + ti*yr£ + \ j y (o° 2 - of)r/ + j a c%r? . (48) 

The invariants l|47|l and (|48|l appear with the independent coefficients W\ and 
w 2 in Q . For the z-component of the field we give first a fuller expression derived 
from a number of low-order invariants 

H z = Ul J z + u 2 T z a +v 1 J z Ol + v % G xy T xyz 

+ Wl (-O zx T£ + O yz Tf) + w 2 2 2 Tf + ... (49) 

Terms in the first line are needed to recover the results by 7 . However, we 
are now only interested in the interplay of T§ octupoles and fields, therefore we 
omit from l|49|l terms not containing 7~f 

H z = Wl (-O zx T£ + O yz rl 3 ) + w 2 0\Tj? . (50) 

If quadrupolar interactions induce any of the quadrupolar moments appearing 
in the above equation, the field will induce T$ octupoles, thus explicitely breaks 
the symmetry of the problem. However, in the absence of such interactions, 
we can turn to the high-field expressions (|37[) - l|39|) which give none of these 
quadrupoles. Therefore, symmetry breaking by octupolar ordering is a possibil- 
ity. 

For fields of other orientation, we need also the following relationships 

H x = w 1 {-O xy T y 8 + O zx T?) - ?f(O 2 + Ol)T£ , (51) 

and 

H y = Wl (-O yz Tf + Q xy T£) + ^(O 2 - Ol)T y P . (52) 
For a field in the (111) direction, 



Hm = Wl [{O xy -O zx )ri + (O yz -0*y)Tf + (O zx -O yz )Tf] (53) 
+w 2 



2 

Neither of the quadrupolar coefficients seen above are field-induced. Therefore, 
if there is no quadrupolar interaction to introduce some of them as order pa- 
rameters, H || (111) fields will allow the same kind of T5 octupolar ordering as 
in the absence of a field (remember, though, that the amplitude of the order 
will be gradually suppressed by the field). 
In a (110) field 

Hno = -f0 2 °(T/-r/) 
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The point to note from (|37Jl - is that though (O zx — O yz ) and 0\ are not 
induced by the field, 0% is, and therefore the octupolar component (%f — T&) 
is also field induced. The remaining octupolar degeneracy arises from the fact 
that (TJ? + T@) and T@ do not couple to the field. Since these are associ- 
ated with different terms in the expansion 

O, ( r x + Ty) and Tz ma y order 

independently. 

Finally, we comment upon the case H |j (11c) 



#11 



tui [(O xy - cO zx ) T£ + (cO yz - O xy ) + (O zx - O yz ) Tf] 



+W 2 



Ol) r% + i (Og - Ol) T/ + cO\T? 



(55) 



Again only (Zjf — T&) is field induced. However, once octupole-octupole inter- 
action gives rise to order, it induces 0\, which in turn induces (T^ + T^), 
thus there is a single phase transition (Fig. right). 



6 Conclusion 

The Tg subspace supports a variety of competing order parameters. The fourfold 
ground state degeneracy can be lifted either in two steps (removing Kramers 
and non-Kramers degeneracies separately) , or in a single phase transition. The 
latter possibility is realized by the ordering of octupoles. We have found 
that the crystal field gives rise to a peculiar single-ion octupolar anisotropy 
which makes the choice of Tf n — T£ + + Tf octupoles preferable as order 
parameter. Though it breaks time reversal invariance, octupolar ordering is 
non-magnetic in the sense of yielding vanishing dipole moments. On the other 
hand, the ordering of T§ octupoles induces quadrupoles as secondary order 
parameter; this feature allows the simultaneous lifting of Kramers and non- 
Kramers degeneracies. 

Our discussion is mainly about a hypothetical Ts-type ferro-octupolar or- 
dering in a lattice of Tg shells. The mean-field results presented in Fig. [2] and 
Fig. |3 are equally valid for ferro-octupolar, and the triple-fc antiferro-octupolar, 
ordering patterns, but our main interest is in magnetic field effects, and our ar- 
guments for the case of non-zero magnetic field are restricted to uniform phases. 

Magnetic octupoles are not time reversal invariant, thus we might have ex- 
pected that spontaneous symmetry breaking due to uniform octupolar ordering 
is necessarily suppressed by magnetic fields. Indeed, for fields of a general ori- 
entation, we find that the degeneracy of different octupolar ground states is 
immediately lifted, and the para-octupolar state prevails at all temperatures. 
However, the analysis of field-induced multipoles shows that for field directions 
lying in certain planes, the field does not induce all the T^-type octupolar mo- 
ments, and therefore sharp octupolar transitions remain possible up to a certain 
critical field. The size of the critical field, and the nature of the transition to 
the high-field dipolar state, depend on the details of field orientation within the 
symmetry-specified planes. 
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